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Abstract. Wc give an asymptotic expression for the number of nonsingular 
integer nxn- matrices with primitive row vectors, determinant k, and Euclidean 
matrix norm less than T, as T — > oo. 

Wc also investigate the density of matrices with primitive rows in the space 
of matrices with determinant k, and determine its asymptotics for large k. 



1. Introduction 

An integer vector v 6 Z" is primitive if it cannot be written as an integer 
multiple m / 1 of some other integer vector w G Z". Let A be an integer n x n- 
matrix with nonzero determinant k and primitive row vectors. We ask how many 
such matrices A there are of Euclidean norm at most T, that is, \\A\\ < T, where 



||A|| := yl>y = s/tv{A t A). Let N' nk {T) be this number (the prime in the 

notation denotes the primitivity of the rows), and let N n ^(T) be the corresponding 
counting function for matrices with not necessarily primitive row vectors. We will 
determine the asymptotic behavior of N' nk {T) for large T, and investigate the 
density D n {k) := limr_ 5 . 00 N' n k (T)/N n ^{T) of matrices with primitive vectors in 
the space of matrices with determinant k ^ 0. 

Let M n ^k be the set of integer n x n- matrices with determinant fc. Then N, lt k(T) = 
\Bt H M„.fc|, where Bt is the (closed) ball of radius T centered at the origin in the 
space M n (R) of real n x n-matrices equipped with the Euclidean norm. Throughout, 
we will assume that n > 2 and k > unless stated otherwise. 

Duke, Rudnick and Sarnak IDRS93^ found that the asymptotic behavior of N U: k 
is given by 

N n , k {T) = c^fcT"^- 1 ) +O e (T n (n-l)-l/(n+l)+ c j ) 

as T — > oo, for a certain constant c n> k and all e > 0, where the error term can be 
improved to 0(T 4 / 3 ) for n = 2. The corresponding case for singular matrices was 
later investigated by Katznelson, who proved in |Kat93j that 

N n , (T) = c nfi T n( - n -V logT + 0(T n( - n -^). 

See the next page for the constants c n ^ and c n ^. 

Let M' k be the set of matrices in M n ^ with primitive row vectors. Then 
N' n k {T) = \Bt n M' n k \. Wigman |Wig05| determined the asymptotic behavior of 
the counting function \Gt H M' n |, where Gt is a ball of radius T in M n (R), under 
a slightly different norm than ours. The results can be transferred to our setting, 
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whereby we have 

Ko(T) = c^oT"^" 1 ) logT + OfT"^- 1 )), n > 4, 

Ko( T ) = C3^ 3(3_1) logT + 0(T 3 ^ log log T), 
^2,o(r) = 4 T 2 ( 2 - 1 )+0(T). 

The case n = 2 above is equivalent to the primitive circle problem, which 
asks how many primitive vectors there are of length at most T in Z 2 given any 
(large) T. 

The main result in our paper is the following asymptotic expression for the 
number of nonsingular matrices with primitive row vectors and fixed determinant. 

Theorem 1. Let k ^ 0. Then 

N' nk {T) = c' nM T n{n -^ +0e ( r »(«-i)-i/(2«)+e- ); 

os T — »• oo for a certain constant c' n k and all e > 0. 

Section [3] is dedicated to the proof of this theorem. 
The constant in Theorem [1] can be written as 

f\YA9)( d - 



,,,n-i 11 Z_,- v a 

I'M d 1 ...d„=|fc| i=l s ]di vy 

for 7^ 0, which may be compared to the constants obtained from [DRS93] . |Kat93| 
and |Wig05| , namely 

C " 
di-<i„=|fc| t=l 

^ n-1 



C( 



n 



C(2) ( " = 2) 

where C is the Riemann zeta function, /i is the Mobius function, and Co and C\ 
are constants defined as follows (these depend on n, but we will always regard n as 
fixed). Let v be the normalized Haar measure on SL„(R). The measure w below is 
obtained by averaging the n(n — l)-dimcnsional volume of E n A u over all classes 
A u := {A G M n (R) : Aw = 0} for nonzero u G K n . In Appendix lAl we give a precise 
definition of w and calculate w(B\). 

Write V n for the volume of the unit ball in R™ and 5„_i for the surface area of 
the (n — l)-dimensional unit sphere in 1™. Then 

K(„-l)Sn-l _ 7T" 2 /2 



Co := w{B 



,'n\ „ ( n(n — 1) 

r - r -i— — i + i 



.. v{Bt n SL„(M)) K(„-i)^«-i Co 
Ci : — lini 



t^oo j<n{n-i) 2C(2)---C(n) C(2) ■ • ■ C(n) ' 

1.1. Density. It will be interesting to compare the growth of N' k to that of N n k- 
We define the density of matrices with primitive rows in the space M nt k to be 

D n Ik) : = lim : — — ^- = . 
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The asymptotics of N n fi and N' n are known from [Kat93| and |Wig05| , and 
taking their ratio, we see that 

for n > 3. We will be interested in the value of D n (k) for large n and large k. The 
limit of D n {k) as k — > oo does not exist, but it does exist for particular sequences 
of k. 

We say that a sequence of integers is totally divisible if its terms are eventually 
divisible by all positive integers smaller than m, for any m. We say that a sequence 
of integers is rough if its terms eventually have no divisors smaller than m (except 
for 1), for any m. An equivalent formulation is that a sequence (k\, k2, . . .) is totally 
divisible if and only if \ki\ p — > as i — > oo for all primes p, and (fci, k%, . . .) is rough 
if and only if \ki\ p — > 1 as i — > oo for all primes p, where \m\ p denotes the p-adic 
norm of m. 

We state our main results about the density D n . We prove these in section |U 

Theorem 2. Let n > 3 be fixed. Then D n is a multiplicative function, and D n (p m ) 
is strictly decreasing as a function of m for any prime p. We have 

D„(0) < D n {k) < D n (l) 

for all k ^ 0,1. Now let fci, &2, . . . be a sequence of integers. Then 

D n {h) -> 1 

if and only if (fci, k2, ■ ■ •) is a rough sequence, and 

if and only if (ki,k2, ■ ■ •) is a totally divisible sequence. Moreover, D n (k) — > 1 
uniformly as n — > oo. 

We prove Theorem [2] for nonzero fc, , but it is interesting that this formulation 
holds for k = also. The case of k — was proved by Wigman |Wig05| , where he 
found that D n (Q) equals l/((n — 1)". We remark that Theorem [5] implies that 

D n (ki) D n {0) 

if and only (fci, ki, ■ ■ ■) is totally divisible, for any fixed n > 3. 

For completeness, let us state what happens in the rather different case n = 2. 

Proposition 3. Let n = 2. Then D n is a multiplicative function, and D n (p m ) is 
strictly decreasing as a function of m for any prime p. We have 

D 2 (ki) -> 

if and only i/lim^oo X^pifc V-P — ^ °°- Moreover, 

D 2 (ki) -> 1 

«/ and onZy j/limj-^oo X^pifej 1/p - ^ 0. T7ie sums are taken over all primes p which 
divide fc,;. 

1.2. Proof outline. Our proof of Theorem Q] uses essentially the same approach 
as |DRS93j . The set M' n k is partitioned into a finite number of orbits ASL n (Z), 
where A £ M n ^ are matrices in Hermite normal form with primitive row vectors. 
We count the matrices in each orbit separately. The number of matrices in each 
orbit scales as a fraction l/fc™ -1 of the number of matrices in SL„(Z). We can 
view SL„(Z) as a lattice in the space SL n (R), and the problem is reduced to a 
lattice point counting problem. The lattice points inside the ball Bt are counted 
by evaluating the normalized Haar measure of Bt n SL„(M). 
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2. Preliminaries 
The Riemann zeta function £ is given by 

oo 1 1 

«-) : =E^=n T 3 1 7^ 

n— 1 p ' 

for Res > 1, where we use the convention that when an index p is used in a sum 
or product, it ranges over the set of primes. 

The Mobius function fi is defined by fi(k) := (—1)'™ if fc is a product of m 
distinct prime factors (that is, k is square-free), and fi(k) := otherwise. We 
note that fj, is a multiplicative function, that is, a function / : N* — > C defined 
on the positive integers such that f(ab) = f(a)f(b) for all coprime a, b. 

We will use the fact that SL n (K) = M n> i has a normalized Haar measure v which 
is bi-invariant (see |Sie45j '). 

2.1. Lattice point counting. Let G be a topological group with a normalized 
Haar measure vq and a lattice subgroup r C G, and let Gt be an increasing 
family of bounded subsets of G for all T > 1. Under certain conditions (see for 
instance |GN10j ). we have 

|G T nr| ~ vq (Gt n G) , 

where we by f(T) ~ g(T) mean that f(T)/g(T) — > 1 as T — > oo. In this paper, we 
are interested in the lattice SL„(Z) inside SL„(K), and the following result will be 
crucial. 

Theorem 4 f [DRS93j . Theorem 1.10). Let Bt be the ball of radius T in the space 
M n (ffi) of real n x n-matrices under the Euclidean norm \\A\\ = y^tr(A t A). Let v 
be the normalized Haar measure of SL„ (M) . Then 

\B T n SL„(Z)| = v(B T n SL n (R)) + o e (r™("- 1 )- 1 /(»+i)+^ 

for all e > 0, and the main term is given by 

\b t n sl„(z)| ~ drf"- 1 ', c 1 = — — — — - — r . 

C(2) • • • C(») r (-) r T ^ + 

In fact, a slightly more general statement is true. We can replace the balls Bt 
in Theorem 0] with balls under any norm on M n (M), and the asymptotics will still 
hold, save for a slighty worse exponent in the error term. 

Theorem 5 ( [G NlOj . Corollary 2.3). Let || • ||' be any norm on the vector space 
M n (M), and let Gt be the ball of radius T in M n (M.) under this norm. Let v be the 
normalized Haar measure of SL„ (R) . Then 

\G T n SL„(z)| = v(G T n SL n (R)) + 0e ( T «( n - 1 )- 1 /( 2 ™)+ £ ) 

for all e > 0. 

We will be interested in the following particular case of Theorem^ Let A € M n ^- 
Then \\X\\' := defines a norm on M n (R), and the ball of radius T in M n (R) 

under the norm || • ||' is A ■ Bt- 

Corollary 6. Let A G M n ^. Then 

\AB T n SL„(Z)| = u{AB T n SL„(R)) + o^^-^-i/CanJ+e) 

for all e > 0, using the notation from Theorem^ 
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3. The number of matrices with primitive rows 

In the present section, we will prove Theorem [1] We begin by noting that the 
divisors of each row in an integer n x n-matrix A are preserved under multiplication 
on the right by any matrix X £ SL n (Z). In particular, if each row of A is primitive, 
then each row of AX is primitive, for any X G SL„(Z). So we get: 

Lemma 7. If A e M' n k then AX G M' n k for all X G SL n (Z). Thus A • SL„(Z) C 
Consequently M' k may be written as a disjoint union of orbits of SL„(Z): 

<,fc= (J ASL n(%), 
AgA 

for properly chosen subsets A of k . In fact, as we will show in the following, 
the number of orbits is finite, and so we may take A to be finite. 
A lower triangular integer matrix 

fc n • • ■ \ 

£, . C 2 1 C 2 2 ' ■ 

: '•■ 

is said to be in (lower) Hermite normal form if < cxx and < c,j < cu for all 
j < i. The following result is well-known. 

Lemma 8 ( Coh93] . Theorem 2.4.3). Assume k > 0. Given an arbitrary matrix 
A G M nt k, the orbit ASh n (Z) contains a unique matrix C in Hermite normal form. 

We may thus write 

m 

M; fe = U^SL n (Z), 

i=l 

where Ax, ... , A m are the unique matrices in Hermite normal form with primitive 
row vectors and determinant k, and m := fe / SL„(Z)|. By counting the number 
of matrices in Hermite normal form with determinant k > 0, we get 

(9) |M n , fe /SL n (Z)|= ^ 

where the sum ranges over all positive integer tuples (d\, . . . , d n ) such that di • ■ ■ d n = k. 
Proposition 10. Let k > 0. Then 

lM; fe /sL n (z)i= ^ IIEm.9) f 

d 1 ...d n =fei=i 9 |d i \y/ 

where the first sum ranges over all positive integer tuples (di,...,d n ) such that 
di ■ ■ ■ d n = k. 

Proof. We want to count those matrices in Hermite normal form which are in M' n k , 
that is, n x n-matrices in Hermite normal form with determinant k and all rows 
primitive. The number of such matrices is 



|M; fc /SL„(Z)| = J2 



d-\---d n — k i—l 



where Vi(d) is the number of primitive vectors (xi, . . . , Xi—±, d) such that < 
Xi, . . . , < d. There is a bijective correspondence between the primitive vectors 
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(xi, . . . ,Xi-i,d) and the vectors y = (yi, . . . ,Vi-i) such that 1 < j/i, . . . , < d 
and gcd(y) is coprime to d. Let d = 1 • • -p" 3 be the prime factorization of d. The 
number of vectors y which are divisible by some set of primes P C , . . . , pj } is 

/ \ i-i 



so by the principle of inclusion/exclusion (see |Sta97] L we have 

i-l 

Vi(d) 



E (-W 

pi,— ,pj] \ 



□ 



We are now ready to derive the asymptotics of N' n k (T) . 

Proof of Theorem [2J Let us write A\ , . . . , A m for all the n x n- matrices in Hermitc 
normal form with determinant k, where m '■= \M' n fc /SL„(Z)|, and let 1 < i < m. 
Then 

\B T nAiSL n (Z)\ = |A l (4ri J B T nSL„(Z))| = | J 4r 1 B T nSL„(Z)|, 
which by Corollary [5] is equal to 

v^Bt n SL„(Z)) + £ (T™(™- 1 )- 1 /( 2 »)+ e ) 

for any e > 0. Since Ai/k 1 / 71 G SL n (R), we get by the invariance of the measure v 
that 

v{A~ x B T n SL n (Z)) = v (A-'Bt n SL n (R))) = 

i/ (V 1 /™^ n ^ SL n (R)) = i/ (S r/fc i/„ n SL n (R)) . 
By Theorem 21 the last expression is equal to 

and thus 

(11) \B T r\Ai SL n (Z)\ = ^r'"" 1 ' +0e ( T n(n-l)-l/(2n)+ E )_ 

Now, 



<, fc (T) = |B T nM; fe | = 

so applying (fTTj) we get 



B T n\jA t SL n (Z) 



J2\B T r\Ai SL„(Z)| 



/"f 

_/V' — 1 jm(n-l) _|_ Q e ^jin(n-l)-l/(2n)+e^ _ 

i=l 

|M;. fc /SL„(Z)| ^L_ T »(«-1) +0e ( T »(«-l)-l/(2n)+e) 5 

and we need only apply Proposition [TU] to get an explicit constant for the main 
term. This concludes the proof. □ 



COUNTING NONSINGULAR MATRICES WITH PRIMITIVE ROW VECTORS 



7 



4. Density of matrices with primitive rows 

Set 

(12) a n (fc):=|M„ jfc /SL„(Z)|= ]T d\---d^\ 

d\ ■ ■ -d n —k 

(13) a' n (k):=\M^ k /SL n (Z)\= £ IIE^Mt) ' 

We would like to calculate the density of matrices with primitive rows in M n ^ 
for k 0, that is, the quantity 

D ( k ) = Km K ' k{T) = <fc = IK. fc /SL»(Z)| = a^M 
n( > T^LN ntk {T) c n M |M„, fc /SL n (Z)| a n (k)' 

We will prove in section T4.ll that a n ,a' n and D n are multiplicative functions, and 
therefore we need only understand their behavior for prime powers k = p m . We will 
now prove a sequence of lemmas which we will finally use in section 14.21 to prove 
Theorem [2] 

Lemma 14. The functions a' n and a n are connected via the identity 



i=0 x 7 

for primes p and m > 0. 

Proof. a n (p m ) counts the number of n x n- matrices in Hcrmitc normal form with 
determinant p m , whereas a' n (p m ) counts the number of such with primitive rows. 
If A is a matrix that a n (p m ) counts which a' n (p m ) does not, then some set of rows, 
indexed by S C [n] := {1, . . . , n} (where \S\ < m), are divisible by p. The number 
of such matrices is a ra (p m ~' ,s '), and thus by the inclusion/exclusion principle, 



4b m )= E (-i) |s| a n r |s| ) = EH) , C «»r')- □ 



SC[n] i=0 

\S\<m 



Lemma 15. For any prime p and m > 1, the following recursion holds: 
a n (p m )=P n - 1 a n (p m - 1 ) + o»-i(p m ), 

or equivalently, 

a »(:P ) = ■ 

Proof. We split the sum 

into two parts, one part where d n is divisible by p, and another part where it is not 
(so that d n = 1)- The terms corresponding to d n = 1 sum to a n -i{p m ). Where d n 
is divisible by p, we can write d n ='■ pe n for some e n . Let := for all i < n. 
Thus, 

E d ?---<" 1= E e?---(e„/p) n - 1 = -^ T M^- 1 ). 

di---d„=p m ei— e„=p m_1 

Adding the two parts gives us a n (p m ) = p n ~ 1 a n (p m ~ 1 ) + a n _i(p m ), from which the 
claim in the lemma follows by rearrangement. □ 
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Lemma 16. Let n and p be fixed, where n > 3 and p is a prime. Then 



A^)^(i-^t)" 



as m — >• oo. 

Proof. We apply the simple upper bound 
a n - 1 (p m )= E 4--<~ 2 < E (P m ) n - 2 = (rn+ir-\p m ) n - 2 

di---d n -i—p m d\ ■ ■■d n -i—p rn 

to the expression for a n (p m ~ 1 ) in Lemma 1151 

anip" 1 - 1 ) = -4zi{a n (p m ) - a n -x(p m )) 

= 4~MP m ) + 0((p m r- 2 (m + l)"- 1 ). 
Repeated application (at most n times) of this formula yields the asymptotics 
an(p m ~ l ) = j-^Ty^n( P m ) + 0((p m ) n - 2 (m + l)"' 1 ) 

for 1 < i < n. 

Now let m — > oo, so that we may assume m to be larger than n. The sum in 
Lemma [T4l then extends up to i = n (because the factors (") vanish for larger i), so 

a' n (j> m ) = j2{-iy( n \ n {p m - 1 ) 

= D" 1 )' (") j-^ryMp m ) + o{{p m ) n - 2 (m + i)- 1 ). 

We divide by a n (p m ) on both sides and use the fact that a n (p m ) > (p m ) n ~ 1 , so 
that 

'n\ 1 /(p m ) n - 2 (m + l) n ~^ 



i=0 

' II 



ij) 171 ) n— ^ 



£(")(^y+^ (m+ir 



% j \p /i j \ p" 1 

V p n - x ) v p m 

As m — > oo, the second term on the right vanishes. □ 

4.1. Proof that D n (p m ) is strictly decreasing. In this section we will prove the 
following proposition. 

Proposition 17. The function D n is multiplicative, and D n (p m ) is strictly de- 
creasing as a function of m for any fixed prime p and dimension n > 2. 

We may rewrite (fT12j) as 

a„ = * • • • * (-) 

where (•)* is the function x M> x % and * denotes the Dirichlct convolution. Similarly, 
we may rewrite (|13J) as 

a^ow-r- 1 )*---*^-) ), 

so by the commutativity and associativity of the Dirichlct convolution we have 

a' n = fj,* n * a n , 
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where fj,* n denotes the convolution of p. with itself n times (so that /i* 1 = /i). Since 
the Dirichlet inverse of \i is the constant function 1, we have also the relation 

a n = 1*™ * a' n . 

As [i and (•)* are multiplicative functions, it follows that a n ,a' n and D n are multi- 
plicative as well. 

Now, we want to show that 

D [ m| _ a 'n(p m ) 

is strictly decreasing as a function of m, for fixed n > 2 and primes p, or equivalently 
that 

(18 ) > <(p m+1 ) 

a n (p m ) a„(p™+ 1 ) 

for all m > 0. 

4.1.1. The case m < n. The inequality (|18p is equivalent to 

<(p m ) > <{P m+1 ) 



(i* n * Ob" 1 ) (i* n *<)(p m+1 ) 

for all m > 0, which is equivalent to 

<{p m ) > <(P m+1 ) 



ET=o l* n (p l W n (P m - 1 ) £™+* l* n ij> l )a' n { P m + 1 - 1 ) ' 
or, after taking the reciprocal of both sides, 

(19) y 1*^)^—1 < y i"V) a "w 

s <(p m ) s <(p m+1 ) 

Now, a matrix A £ M' n pm in Hermite normal form with primitive rows can be 
mapped uniquely to a matrix A' £ M' n pm+ i by multiplying the largest diagonal 
element by p, where we break ties by always choosing the diagonal element on the 
later row. This shows that m n- a' n (p m ) is a non-decreasing function, and thus 

all factors ^gfe^n - and a „/^ p m+i^ be in the interval (0, 1], and since also 1*™ is a 
positive function, the last inequality Q19p holds if 



i=0 



i*™(j/) < r n {p m+1 ). 



Now, \* n (p l ) is the number of ways of writing i as a sum of n non- negative integers, 

n-l+i 
ra-1 



n — 1 + i\ ^ ( n + m 



and it is well-known (see |Sta97j ) that this is equal to ("^J 1 ^ 4 ) for all i, so the last 
inequality is equivalent to 

n l 

y, 

^— ' » n — 1 J \n— 1 

A well-known combinatorial identity (see exercise 2.1 in |Sta97| ) states that the 
left side above is equal to ( n ~ i n m ) ■ By the unimodality and symmetry of binomial 
coefficients, we have ("+"') < (™+™) if and only if (n + m)/2 < n - 1/2, which is 
equivalent to the inequality m < n — 1. We have therefore proven (|18p and thus 
Proposition IT71 for m < n— 1. 
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It remains to prove (|18[) for m > n. We begin by making the following obser- 
vation. The inequality (fT5|) is equivalent to (fl^j) . and the inequality (fT!))) holds for 
m > n if 

<(p m "') < <(p m+1 - 1 ) 
<(p m ) ~ <(P m+1 ) ' 
for all i < m. We can rearrange this inequality as 

<{p m+1 ) < <(p m+1 - i ) 

which states that 

<(p m+1 ) 
<(p m ) 

is a non-increasing function of m > n, for fixed 77 > 2 and p prime. We will therefore 
be done if we can prove that 

(20) a' n (p m )a' n (p m+2 ) < a' n {p^)a' n {p m + l ) 
for all m > n. 

4.1.2. The case m > n. We will now prove (|20p . Accordingly we will assume 777, > n. 

We begin by noting that a n (p m ) can be written as the Gaussian binomial coeffi- 
cient (see |Sta97j ) 

/m + n-l\ (p^ 1 -!)■■■ - 1) 

(21) ^ )= { n-l )= (p-D-.fcn-i-U 

which may be proved by simply observing that ( m ^" 1 ~ 1 ) p = P™" 1 ( (m_ ^[ n-1 ) + 
/m+(ri-i)-i\ sa ti s fi es t;h e recursion formula for a n {p m ) given in Lemma [T51 with 

the same initial values, and thus must coincide with a n (p m ). 

Now, the numerator a n (p m ) ■ (p — 1) • • ■ (p ri_1 — 1) of (|21j) equals 

n-l 

(22) (p m+1 - 1) • • ■ (p"+»-i - l) = 5^(-l) n - 1 - i j» rai Qi(p) 

i=0 

where for each z we have defined the polynomial 

Qi(p):= E P C1+ - +C % 

l<ci<— <c 4 <n— 1 

which depends on n but not on m. Thus, using the formula for a! n {jp m ) from Lemma 
[Til we get 

•(p-l)...(p»- 1 -l) = 

N n-l 

7l N 



E(- X n ; )E(- 1 ) n ~ 1 " i f (m " i)i< 2«(p) 

j=0 i=0 

n— 1 n / \ 

i=0 j=0 ^ ' 

n-l 

£(-i) n " 1-i Qi(p)( 1 -p" i )V n4 . 



i=0 
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where we used the fact that for m > n, the sum over j extends up to j = n. 
Inserting the expression above into (|2T))) . we find it remains to prove that 

n—l n — l 

i=0 j=0 
n—l n—l 

i=0 j=0 

This can be rewritten as 

n-l 

E (-i) i+J 'Qi(p)Q»(i -p~T(i - p -Tp mi+{m+2)3 < 

i,j=0 
n-l 

E (-i) i+i Qi(p)^(p)(i-p" i )"(i-f" i )V m+1)(i+i) 

or equivalcntly 

n-l 

< E (-l^'QiCpJQjCpJCl -p"T(l -p-J) n (p( m+1 )( i+ J) _ p ™+(™+2b). 

Since p( m + 1 )('+i) - p «»+(™+2)j = _ p m{i+j) p 2j^ the last i neqU ality is 

equivalent to 

n-l 

o < E (-i)^Q l (p)g J (p)(i -p^)V" (?;+J) (p' i+J -p 2j )- 

We observe that the diagonal terms i = j vanish, as well as the terms with i = 
and the terms with j = 0. Pairing terms (i,j) and opposite the diagonal, and 
observing that 2p t+: > — p 2% — p 2 - 7 = — (p l — p 3 ) 2 , we get the equivalent inequality 

(23) o< E (-i) 1+l+3 Q t (p)QMa-p'Ta-p- 3 ) n p m{l+j) (p l -p 3 ) 2 - 

l<i<j<n-l 

for all m > n och all primes p. The sum is empty for n = 2, so we may assume 
that n > 3. We note that each term in the sum (|23|) is positive if i + j is odd, 
and negative if i + j is even, and thus each negative term satisfies i + 1 < j. We 
will show that the sum is non-negative by showing that each negative term with 
position is no larger in absolute value than the positive term with position 

(i + 1, j). That is, 

Q 1 {p)QM^ -p"T(i -p- ] ) n p m{l+1 \p l -p 3 ) 2 < 
Q l+ i{p)QM^-P^ l Y\^-P^ n P m(l+1+j) {p l+1 -P 3 ) 2 , 

which is equivalent to 



Q t (p)(i~ P -T(p t -p r ) 2 < 
q, +1 (p)(i-p- 1 )T(p' +1 -p') 2 
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We prove the last inequality by comparing the left and right-hand sides factor by 
factor. We have i + l<j<n — 1, and thus i + 1 < n — 1, and therefore 



1<C1<-<C;<?1-1 

l<ci<-<c i+ i<?i-l \r r r / 

l<Ci<---<Ci_|_i<?l— 1 

since both sums are nonempty and 1/p-h l/p 2 H < 1 for all p > 2. We also have 

(i-p-T ^(i-f^- 1 )"- 

Finally the factor (p> — p 1 ) 2 is smaller than the factor p m {jp — p l+1 ) 2 , which is 
demonstrated by the following string of implications 

p(p - I) 2 > 1 =^ p 3 (l - l/p) 2 > 1 => 

p m {\ - l/p) 2 > 1 p m {l-p l+1 /p>) 2 > 1 ^ 

p m {p> -p l+1 ) 2 > {p>) 2 =^ p m (p> -p l+1 f > (p> -p 1 ) 2 , 

which hold for all m > n > 3, j > i + 1 and p > 2. Thus (f2"0]) is true for all m > n. 
This concludes the proof of Proposition (fTTl) . □ 

4.2. Asymptotics of the density function. In this section we prove Theorcm[5] 
and thus derive the asymptotics of D n (k). Fix n > 3. For any nonzero integer ki, 
write fc.j = YlpP 171 "^ as a product of prime powers, where all but finitely many of 
the exponents m p (i) are zero. Then since D n is multiplicative, we have 

D n (h) =l[D n (p m ^). 

p 

Now, by Lemma [T51 and Proposition [T71 we get 

i>n c „(p^>)>n(i-^)"=^>o, 

so it follows that Y\ p D n {p mp W) is uniformly convergent with respect to i, and 
therefore 



(24) lim TT D„ (p m ^) = TT lim D n (p 



P P 



Let (fci, k2, . . .) be a sequence of nonzero integers. It now follows from (|24[) . Propo- 
sition [T7] and the fact that D n (l) = 1, that 

D n {h) -> 1 

if and only if m p (i) — >■ as i — > oo for all p, that is, if and only if (ki, k,2, ■ ■ • ) is a 
rough sequence. Likewise it follows, using Lemma (|16p . that 

D n (ki) -> — ; — - — — 
»V ^ ^ n _ ^„ 

if and only if m p (i) — > oo for all p, that is, if and only if (fci,/c2, ■ ■ ■) is a totally 
divisible sequence. Since D n (0) = 1/C( n — we may allow the elements of the 
sequence (fci, fe, • • ■) to also assume the value 0. 

Finally, it follows that D n (k) — > 1 as n — > oo uniformly with respect to k since 
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as ri — y oo because £(n — 1) = 1 + 0(2 n ) for n > 3. We have thus proved all parts 
of Theorem [2 □ 
We conclude this section by proving Proposition^ which tells us the asymptotics 
of D 2 (k) for n = 2. 

Proof of Proposition^ If m = 0, we have D2(p m ) = 1. Assume m > 0. The 
2 x 2-matrices in Hermite normal form with determinant p m and primitive rows 

are of the form \ ® m ] where < x < p m ,p \ x. Thus a' 2 (p m ) = p' n (l - l/p). 



Moreover, 



x p 



m r,m+l 1 1 1 /„m+l 



d\d2—p m i-\-j—m i—0 

so D 2 (p m ) = (1 - l/p) 2 /(l - l/p m+1 ). Therefore 

1\ 2 1 
1__ <£ 2 (p™) < 1--. 

P/ P 
Since D 2 is multiplicative, we get 

2 



< Daw < n ( 1 

p\k 



The left and right sides both tend to if and only if lim^oo J2 P \kt VP — ^ °°> an d 
they both converge to 1 if and only if lim^oo J2 P \kt l/p — >• 0. □ 

Appendix A. Calculation of a measure 
In |Kat93j the asymptotics 

N, h0 (T) = ^—±w{B)T n ^ logT + 
CO) 

are given, where -B is the unit ball in M„(R). The measure iu on M n (K) is de- 
fined in |Kat93j as follows. Let A u := {A e M„(R) : = 0} be the space of 
matrices annihilating the nonzero vector u G M" \ {0}. Wc define for (Lebcsguc 
measurable) subsets E C M n (M) the measure w u (E) := vo\(E n A u ) where vol is 
the standard n(n — l)-dimensional volume on A u , and define the measure w(E) := 
(1/2) Jg„_i w u (E) dv[u), where f is the standard Euclidean surface measure on the 
(n — l)-dimensional sphere S Il_1 . 

We shall now calculate w(B). The set B D A u is the unit ball in the n(n — 1)- 
dimcnsional vector space A u . Its volume does not depend on u ^ 0, and if u = 
(0, . . . , 0, 1), then B n A u is the unit ball in M™^" 1 ), when identifying M„(R) with 
M" 2 . Denote by V^n-i) the volume of the unit ball in R"(" -1 ). Thus w u (B) 
1 independently of u 7^ 0, and 

w(B) = K (n _i)- / _ dv{u) = , 

where SVi-i is the surface area of the sphere S™ _1 . The volume and surface area of 
the unit ball is well known, and we may explicitly calculate 



Co := w(B) = 



-n-" 2 /2 



2, 
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Recalling from Theorem 2] the expression for C\, we get the following relation. 

r 1 1 r 

1 C(2)-C(") r pj r f "("-1) , ^ C(2)---C(n) °- 
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